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PACS 72 . 25 . Rb - Spin relaxation and scattering 

PACS 85.75. -d - Magnetoelectronics; spintronics: devices exploiting spin polarized transport or 
integrated magnetic fields 

Abstract. - We consider the magnetic response of a two-dimensional electron gas (2DEG) with 
a spin-orbit interaction to a long-wave-length electromagnetic excitation. We observe that the 
transverse electric field creates spin polarization perpendicular to the 2DEG plane. The effect is 
more prominent in clean systems with resolved spin-orbit-split subbands, and reaches maximum 
when the frequency of the wave matches the subband splitting at the Fermi momentum. The 
relation of this effect to the spin-Hall effect is discussed. 
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Spintronics is a fast developing field of research. While 
the spin-orbit effects in semiconductors have been dis- 
cussed for a long time [1-3], recently discovered spin-Hall 
effect [4,5] has generated a lot of interest in applying spin- 
orbit related effects to spintronics [6]. The promise of 
the spin-Hall effect is in the possibility of generation of 
non-equilibrium spin polarization by means of a DC elec- 
tric field. However, after some discussions it was con- 
cluded that in two-dimensional electron gases with the 
Rashba and the Dressclhaus spin-orbit interactions, the 
spin-Hall effect vanishes for constant and homogeneous 
electric field [7-10]. At finite frequencies the spin-Hall ef- 
fect is non-zero [8]. 

While for the homogeneous spin-Hall effect (uniform ap- 
plied electric field) the out-of-plane spin polarization is ex- 
pected to accumulate only at the edges of the sample [8] , 
there is an alternative possibility that we propose to ex- 
plore here: To create an out-of-plane inhomogeneous spin 
density in the bulk in response to a spatially-modulated 
field. Such bulk accumulation is free from the uncertain- 
ties associated with the charge and spin transport near the 
sample boundaries, and thus may provide an unambiguous 
method to detect the spin-Hall effect. 

A possible realization is shown in Fig. 1, and cor- 
responds to a transverse electric field E x changing in 
the y-direction. In the presence of the spin-Hall ef- 
fect, the spin currents are expected to flow towards each 
other, with the spin density accumulating in the mid- 
dle. From the Maxwell equation, xotE = —(l/c)B, a 




Fig. 1: The "naive" motivation: Spin accumulation due to 
converging spin currents. 



transverse electric field can only exist together with the 
time-dependent magnetic field. To create the situation 
shown in Fig. 1 we therefore consider a linearly polar- 
ized in-plane microwave field A = e x Ao exp(iq y y — iQt); 
E = E x e x = e x (A / c)itt eyLp(iq y y - iCtt); B = B z e z = 
—e z Aoiq y exp(iq y y — iClt); and fl = c\q\ = cq y . There 
is naturally a spin response in the z-direction due to 
the B z magnetic field (the Zeeman effect; this is simi- 
lar to the response to an in-plane field considered, e.g., in 
Rcfs. [11, 12]). However, in addition to this pure spin ef- 
fect, we also find an orbital contribution, g° rbltal rotE, 
which can be attributed to the spin-Hall effect. (Since 
rotE ~ £IB Z this response can be conveniently interpreted 
as an orbital contribution to the spin susceptibility.) The 
response of the in-plane polarization to the transverse elec- 
tric field at arbitrary values of q was considered for the 
clean case in Ref. [13]. 

We consider the Rashba model of a 2DEG with the 
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Hamiltonian 

(P--A) 2 f e \ 1 

H = — — \-anri[P A ) - -^g^Bcr + V disoidci . 

2m* V c / 2 

(1) 

where r\ = z x <x = (—a y ,a x ,0) and m* is the effec- 
tive electron band mass. The model (with A = and 
Vdisordcr = 0) has two energy bands. The eigenfunctions 
of the two bands j3 = ± are given by 



1 



t±>p ) = ^(lt)T^U)) , 



(2) 



where tan< 



= Py/Px and the energies are given by 

e± (P) = ± A p> wliere A p = OrIpI- 

One should distinguish between the clean and dirty 

regimes, a) In the clean regime (weak disorder or strong 
SOI) Ap e a R p f > r -1 , where r is the elastic mean free 
time. The two spin-orbit split bands are, thus, well de- 
fined, b) In the dirty regime (strong disorder or weak SOI) 
Af <C t -1 . The band splitting is not resolved. This is the 
regime where the Dyakonov-Perel relaxation applies [1]. 
One can also define a "super clean" regime 1/r < m*a R , 
which, however, appears to be irrelevant experimentally 
at this time and will not be considered here. 

In this work we focus primarily on the clean regime, 
although our diagrammatic derivation holds also in the 
dirty regime. Depending on the frequency of the applied 
field one can then probe the resonant (Q m 2Ap) and off- 
resonant — 2Ap\ S> 1/r, e.g. f2 <C 2Ap) responses 
of the spin system. Before rigorously deriving our results 
diagrammatically, we provide qualitative derivation of spin 
responses in these regimes. 

In the low-frequency (off-resonant) limit, we consider 
slow (both in time and space) components of the single- 
electron density matrix //9 1 ./3 2 (Pi,P2)- Here Pi and 
are the band indexes while p 1 and p 2 are the electron 
momenta. Introducing p = (p 1 + p 2 )/2 and q = p 1 — p 2 
we note that slowness in space means that \q\ is small. 
When the spin-orbit splitting, A p , is large, only the band- 
diagonal elements fp{p\,P2) = ff3,p{PiiP2) can be slow. 
Indeed the band-off-diagonal elements either oscillate in 
time with frequency ~ A p or oscillate in space with |q| ~ 
Ap/vp. Consider the density matrix of band f3 in the 
coordinate representation: 

ffs(r,r') = 

MPi,P 2 )e^ r - ip ^ \u , Pl )(u^\ .(3) 

P1P2 

We calculate the z-componcnt of the spin density (cf. 
Ref. [14]): 

= ^Tr^/^^r)] 

~ \ Yl fp<i>0 eiqr ( u P>p\ u * q ' V P \ u P>p) 

^ npq 



(4) 



f},pq 



Note that the contributions of both bands add up. For 
q = q y e y and using V Ph p = cos (f> p /p we obtain 
s z (q) oc iq y J2pp fp,q,@ C0S 4>p/P- From the usual Boltz- 
mann equation we know that under the influence of the 
electric field E = E x e x the correction to the distribution 
function satisfies J2p fp.q.p cos 4>p veE x VpT, where v 
is the density of states and r is the momentum relaxation 
time. Thus we finally obtain s z {q) oc iq y eE x r (we have 
used pp = m*VF, v oc m* and the fact that everything 
happens near the Fermi energy). 

When the excitation frequency approximately matches 
the spin subband splitting (the resonant case) , the genera- 
tion of the transverse spin polarization can be understood 
based on the Bloch equations, familiar from the theory of 
the nuclear and electron spin resonance [15]. The spin- 
orbit coupling can be represented as a momentum depen- 
dent magnetic field 60 acting on the electrons with the 
Hamiltonian Hso = b$cr, where 60 = ^p{Py/Pi—Px/p)- 
The spin-dependent part of the interaction of the EM 
wave with electrons can be recast in the following form, 
Hint = bier where 61 is a small oscillating magnetic field 
directed along the y-axis, bi = aR(e/c)A e y e tqyy ~ tnt . 
The single-electron states |u-,p) and \u+ iP +q) that are 
connected by the AC field b\ do not have the same mo- 
mentum, and thus their energy difference is given by 
5E p = e + (p + q) - e~(p) ~ 2A p +p y q y /m*. For each 
such pair of states, provided that, without the driving, 
\u- y p) is occupied while is not, the usual Bloch 

dynamics of the electron spin takes place. That is U-,p) 
plays the role of the state ||) along the direction of 60 , 
while jii+^+q) plays the role of Only the part of b\ 
perpendicular to 60 is important, |&i.j_| = |bi|sin^> p . In 
the linear response regime the driving produces the out- 
of-plane polarization 
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b\ I sin <j) p 



[i/T 2 + (fl - SEp)] 
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where T 2 is the dephasing time. T 2 ~ r in the clean 
limit. The total spin polarization is given by s z = 
L^ 2 ^2 p (o- Zi p)np(e~)[l — np(ep +q )}. The summation is 
effectively restricted to the ring defined by pp + a^m* > 
p > pp — an,m*. (Note that in the limit q = the spin 
polarization vanishes, since (cr Zl p) oc sin0 p , i.e. it has 
only a p-componcnt. This actually means that that there 
is a resonant spin current j* flowing in the system [8].) 
Restricting ourselves to the experimentally relevant "not- 
super-clcan" regime 1/r > a^m* , all states in the ring 
pp 4- aRm* > p > pp — a^m* are "in resonance" si- 
multaneously (cf. (5)) when fl ~ 2Ap. Assuming also 
vp\q\ < 1/r we expand the denominator of (5) and after 
summation obtain 



Rc 



ie- tilt (e/c)A qA 



[i/T 2 + (0 - 2Ap)}< 



(6) 
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Notice the strongly resonant (Fano shape) nature of this 
contribution. A very similar result also obtains for the case 
of the Drcsselhaus spin-orbit interaction. Interestingly, by 
noting that Aq/c = E x /(iCl) and replacing Cl — > ijr we 
can approximately recover the low-frequency spin accumu- 
lation from the resonant result. This indicates the close 
connection between these two effects. 

A resonant generation of out-of-plane spin polarization 
was also recently proposed in Ref [16]. In that work, the 
polarization appeared due to explicit symmetry breaking 
by a strong constant magnetic field applied in the 2DEG 
plane, while in our case the effect is due to the finite value 
of q. 

Diagrammatic derivation. We now derive these results 
in a rigorous manner. We follow the route of Refs. [8,17] 
and derive, first, the kinetic (diffusion) equations govern- 
ing the dynamics of the charge and spin densities. We 
use the Keldysh technique. For introduction see Ref. [18]. 
Although such derivations are standard [18], for complete- 
ness, we provide here the main elements. 

We consider only s-wave disorder scattering, that is 
Vdi Sor dor = y]i- uS(r — rfc) where are random locations 
with the average density 7ii mp . We employ the linear re- 
sponse: H = H + Hi where 



Hn = 



2m* 



anVP + disorder 



(7) 



and 



Hi 



e 1 

Y { v > A }+ + e &- -j9mBa- 

-£{(£ +a *»»M + + * 



(8) 



To make the calculation gauge invariant we have added 
a scalar potential perturbation, <fi. We neglect the A 2 
contribution and introduce, for convenience, a 4-potcntial 
$ = e<M) - (l/2)gfi B B(T. 

In order to account properly for the vertex correc- 
tions [7,9] we perform the linear response expansion start- 
ing from the Dyson equation which reads (id/dt—H)*G = 
1 + S * G , where * stands for convolution. The self-energy 
is given by the self-consistent Born approximation. 



S(ri,r 2 ,ti,t 2 ) = mt 2 <5(ri - r 2 )G{r 1 , r 2 , h, t 2 



(9) 



Using H = H + H 1 ,G = G + G 1 , and S = S + Si we 
obtain 

Gi = G * (#i + £i) * G • (10) 

The inclusion of Si is equivalent to accounting for the 
vertex corrections. This can be seen from Fig. 2. 

The zeroth order in A and Green's functions and 
self energy reflect the standard disorder broadening. One 
obtains S^ = —i/(2r), where r^ 1 ~ 2nni mp u 2 v and 
v = m* /{2irh 2 ) is the density of states (per spin) with- 
out the spin-orbit coupling. The zeroth order (in A and 
$) disorder-broadened Green's functions read 



1 1 r\p 



G n 



lrjp 

2 p 



G 



R- 
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Fig. 2: Diagrammatic representation of Eq. 10. The dashed 
lines represents the sum of impurity "rainbow" diagrams (Dif- 
fuson) with minimum one impurity line. 



(lu - e ± (p) +i/(2r)) . In equilib- 



where Gq ± (p,lu) 

rium Gq = h{oj) (G R - G$) and E{f = h(cu) (S£ - S5 4 ), 
where h(uj) = tanh - ~j? F . 

Our goal is to find Gi with the self-consistency condition 



Si(q,w,n) 



d 2 p 

(2tt)2 



Gi(p,q,u,Cl) 



(12) 



Having Gi one can calculate any single-particle quantity, 
i.e., density or current. Thus this approach is equivalent to 
the Kubo formula with the external vertex (the responding 
observable) being chosen later. 

As usual in the linear response theory the Keldysh com- 
ponent Gi splits into two parts, Gf — Gf" ' 7 + G 7 ^' 77 . 
Accordingly does the self-energy, Sf = sf"' 7 + sf"' 77 . 
The first part, G 7 ^ 7 , corresponds to the R-A combina- 
tions in the Kubo formula, while G 7 ^' 77 stands for the 
R-R and A-A combinations [18]. The spin-charge density 
matrix is given by p = § + ser = f^J [-i G<] = 



I / f I £fr(G? ~ G? + Gf) = p'+ P 77 , where p 1 



2 J 2w J (2tt) 2 

-|J" / (CyjCf' 7 and P U represents the rest. Here 
n(q, Cl) is the charge density while s(q, Cl) is the spin den- 
sity. We obtain 

Gf J {p, q, w, fl) = G R {p + q/2, lu + Cl/2) • 

• {h(tu + Cl/2) -h(tu- Cl/2)} H 1 + ^ J {q,LU,Cl) ■ 

■ G5 4 (p- q/2, a; -0/2) . (13) 

and Gf' 77 = Gf h(w - Cl/2) - h(tu + Cl/2) Gf, where 
the correction to the retarded Green's function is 
given by G R (p,q,co,Cl) = G R {p + q/2,u + Cl/2) ■ 
[#i(p, q, w, CI) + Sf (p, q, u, Cl)] ■ G R (p - q/2, to - Cl/2). 
The dot product is used here to indicate the matrix mul- 
tiplication. 

Our calculations below apply both in the clean (but not 
"super-clean") regime, m*a 2 l < r _1 < Ap, and in the 
dirty regime, r _1 > A^. In the "super-clean" regime, on 
the other hand, the integration over lu in the window of 
width Cl around the Fermi energy (this window is cut by 
the function h(ui+Cl/2) — h(co— Cl/2) in (13)) is not smooth 
and more care has to be taken. 

Integrating Eq. (13) over p and lu we obtain 



= iciui 



e{vA}^ 



2c 



(14) 
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The functional / is defined as 
I[X{P)] = 



1 



d 2 p 



j m"T J (2ir) 2 
G£(p-q/2,E F -n/2) 



G^(p + q/2,E F + n/2)-X(p) 



(15) 



while I is a 4 x 4 matrix defined by its action on the 
4- vector p as Ip = I[p] (we just use the fact that p is inde- 
pendent of p to represent the functional I[p] as a product 
of a 4 x 4 matrix / and a vector p). 

The second contribution to the density, p 11 is given by 



p" = . 
Thus the total density follows from 



(p + v<&) = iVtvI 



e{vA}_ 
2c~ 



- $ 



(16) 



(17) 



We rewrite this equation to have all the terms associated 
with the perturbing fields on the RHS: 



-P 



i^lvl 



e{vA} A 



2c 



[(1 - iCtr)! - 1] 



z/<l . (18) 



We allow for arbitrary external frequency f2, including 



n > 



However, we limit ourselves to the experimen- 



tally relevant regime vf|<z| <C t -1 . (Recently the spin and 
charge response functions to the longitudinal fields were 
calculated for arbitrary values of q in Ref. [19].) Expand- 
ing the RHS of Eq. (18) in powers of q up to the terms 
linear in q we obtain 



(1-^) 



-v(eE ■ rj) 



a R b 2 



(iq x eE) ■ a 

b 2 



2a{a 2 
b 2 (a 2 -b 2 ) 
8na{a 2 + b 2 y 



b 2 ) 



4r(a 2 + b 2 ) 



ugp B (B 



ByOy 



(19) 



where E = iftA — iqcj) (we also assumed divE = 0). We 
have introduced a = 1 — iilr and b = 2 Apr. Note that the 
magnetic terms (oc B) are already of first order in q. The 
term in the RHS of (19) proportional to iq x E — rotE 
will lead to our central result. The q expansion of the 
inverse "diffuson" in the LHS of Eq. (18) reads 



1 - I 

T 
1 

r 
1 

+ - 

T 



1 

1 - - 

a 



n 
2 



1 



2a 2 



where is the diffusion tensor. In the dirty limit, b <C 
1, we obtain D^ n w {D / a?)8 nm 8p 1 , where D = v f t/2 is 
the diffusion coefficient. In the clean limit, b ^> 1, only 
the diffusion coefficient in the charge channel is important 
and we obtain D™ m w {D/a 3 )S nm . Eqs. (19) and (20) 
generalize the results of Rcfs. [8,17] for arbitrary frequency 
Q and field polarization. Expanding these equations in 
the dirty limit, b <C 1, and for fir < 1 we obtain the 
equations of Ref. [17]. The two last terms of Eq. (20) taken 
in the limit Qt <C 1 reproduce the spin-charge and spin- 
spin couplings of Ref. [8] . More importantly, we introduce 
a new driving term in these equations, i.e., the second term 
in the RHS of Eq. (19) proportional to rotE. This term 
was not appreciated in Refs. [8,17] where only longitudinal 
driving fields were considered. 

We will now use these equations to obtain the spin den- 
sity. We calculate contributions of zcroth and of first 
orders in q, s = Sq + s±. The first order contribu- 
tion further splits into the orbital and the Zeeman parts, 

„ ...orbital _j_ ,,Zccman 

* 1 — &i ' ^1 

Zeroth order in q. The zeroth order in q corresponds to 
a homogeneous field. In this case there is no difference be- 
tween the transverse and the longitudinal cases. Thus we 
should reproduce known results. For the spin polarization 
we obtain (cf. Ref. [20]) 

b 2 , s 

s = a K eE x rve y — — — . (21) 

b 2 (2a - 1) + 2a 2 (a - 1) 

For fir <C 1 in the case of strong spin-orbit (6^1) and for 
ft <C A 2 f t for the case of weak spin-orbit (b <C 1) this gives 
Sq(CI — > 0) ~ a R eE x rve y . In other words this result is 
obtained when Q,t s <C 1, where r s = max[r, l/( A f t)] is 
the spin relaxation time. This is the well known result [2,3] 
meaning that there is an in-plane spin polarization per- 
pendicular to the applied electric field. In the clean limit, 
Apr 3> 1, at fi « 2Ap the polarization (21) shows reso- 
nance discussed in Ref. [12]. 

First (linear) order in q. After some algebra we obtain 

^orbital 



a R b 2 



Aa 2 {a 2 + b 2 ) 



&2 

{*7,*<?/S}- 



2a(a 2 + b 2 

cr /3- D m«( i 9rO(^m) Tl 

iabvp 



SyUy) 

pa 1 



2(a 2 + b 2 ) 



2 iVMp]- 



(20) 



ieE x 



4tt 



(b 2 + a 2 - a) 



b 

2a 



To- 



ab 



vpbB z 



b(l - a) 



(b 2 



a) 



b 

27i 



-b 2 
ab 



U 2 



(22) 



where Tn 



2(q-l)(j) +q~) 



= b 2 {2 a-l)+2a>(a-l) - W ° haVC UScd % E * = 

-£lB z /c and ^ = (^f-) ufj, s , where p B = e/(2m e c) 
and m e is the bare electron mass. Thus we obtain the 
out-of-plane spin polarization. This is the main result of 
this work. It follows from the second term in the RHS of 
Eq. (19). For Q.t s < 1 

ieE x Tq y m e iVLr 

a = e * — V ^B Z — 

8tt m* 2 

The out-of-plane polarization is much smaller than the in- 
plane one given by Eq. (21). Indeed we observe that for 



^orbital 



(n) 



(23) 
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Fig. 3: Real part of the total spin susceptibility \ = s z /B z for 
2 Aft = 5. Parameters assumed as in GaAs: m e /m* « 15, 
g « -0.44. 



Or s <C 1 their ratio is given by |s° 



orbital I 



*0 



\q\ho, 



where lso = h/(m*ctTi) is the spin-orbit length. The fact 
that the result is purely imaginary means that the spin 
wave is phase shifted by tt/2 relative to the applied mi- 
crowave. This is qualitatively consistent with the naive 
picture of the z-polarized spin currents proportional to the 
electric field (see Fig. 1). However, for quantitative analy- 
sis one would need to use the correct continuity equation 
—iils z + iq y jy = t z , where t z is the torque density [17]. 

Near the resonance, Vt ~ 2A^, in the clean limit, 
Apr 3> 1, we obtain 



^orbital 



m* 



A 2 

F 



(f7-2A F + ^:)(fi-2A i . 



3* ^ 



(24) 



Thus we obtain a double pole like in Eq. (6). 

Response to the Zeeman term. The Zeeman response 
follows from the last term in Eq. (19) and is given by 



„ Zeeman 



9 



2 b 2 - (1 - a)a 



(25) 



In the limit of zero spin-orbit coupling, i.e., b — 0, but 
finite £1 we observe that g Zeeinan vanishes. The physical 
meaning of this result is that without the spin-orbit cou- 
pling no spin flips are possible and, thus, the spin den- 
sity cannot react to a dynamical perturbation. The other 
limiting procedure (first fl — > 0, then 6^0) gives the 
thermodynamic Pauli susceptibility. 

Discussion. In the clean limit Apr 3> 1 the orbital out- 
of plane spin polarization (22) dominates over s ^ eeman near 
resonance fl ~ 2Ap (although s^ eeman is also resonant 
there). Also in the clean limit the low frequency out- 

" > 



•bitall 



of-phase orbital response is dominant, i.e., Im[s° 
Im[si| eman ] for Or < 1. In the opposite regime, A f t < 1, 
the ratio of absolute values of the orbital and the Zeeman 



Fig. 4: Imaginary part of the total spin susceptibility \ = 
Sz/Bz for 2Aft = 5. Inset: Blow-up of the imaginary part 
at Q.T < 1. Parameters assumed as in GaAs: m e /m* « 15, 
g « -0.44. 



contributions becomes independent of Ap, 

fir 



^orbital I 



1 me 
9 m* 



Vi + n 2 r 2 



(26) 



If the prefactor (l/g)(m e /m*) is large (as it happens, 
e.g., in GaAs) the orbital contribution dominates at high 
enough frequencies. 

It is convenient to measure the susceptibility in units of 
v[ib- In Figs. 3 and 4 the resonant behavior at high fre- 
quency is shown. For plotting we choose the clean regime, 
2Ai?r = 5. The orbital Fano-like contribution is clearly 
pronounced. The linear in frequency dependence of the 
imaginary part at Qt < 1 (shown in the inset in Fig. 4) 
corresponds to Eq. (23). 

For numerical estimates we take the parameters of the 
2DEG in GaAs/AlGaAs heterostructures from Ref. [21]. 
With m* - 0.067TO e ~ 0.7 • l(T 31 kg we obtain the den- 



m*/(2irh z ) 



10 36 J 



With 



23 J/T we obtain ufx B ~ 10 1 
With the mobility /Lt = er/m* 



25m 2 /Vs 



(i 



sity of states v 
Mb - 10" 
10 5 G- X cm 

and the electron density n Q 
t ~ 10 _11 s, vp ~ ^iio/{m*v) ~ 3 • 10 5 m/s, and the 
mean free path I = rvp ~ 3/im. The reported value for 
the spin-orbit length lso ~ 4/iTO. Thus we are on the bor- 
der between the Dyakonov-Perel (dirty) regime and the 
strong spin-orbit regime so that Aft = l/lso ~ 1- 

The oscillating transverse electric field needed to gen- 
erate the bulk spin accumulation described here can be 
created inside a microwave wave guide or a strip-line res- 
onator. Assuming for the estimate a microwave field with 
B z ~ 1 Gauss, the spin polarization on resonance is about 
2 • 10 7 cm~ 2 . Such polarization should be measurable with 
the magneto-optical spectroscopic techniques [22,23]. Al- 
ternatively, the transverse (circulating) electric field could 
be generated with an AC magnet. 
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